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Introductwnjo  ^  ^  0NR-sponsored  research  was  to  obtain  the  phonon  spectra  and  elastic 
constants  for  crystals  of  the  explosives  HMX,  RDX  and  CL20  by  vibrational  spectroscopy.  The 
spectroscopies  used  were  Raman  and  Brillouin  scattering.  Another  objective  was  to  measure  the 
modal  strain  Griineisen  parameters  of  these  systems  using  piezomodulated 

these  data  are  crucial  to  the  development  and  testing  of  current  models  of  detonation  that  are 
concerned  with  the  initial  transfer  of  mechanical  energy  mto  molecular  degrees  of  freedom.  The 
intent  was  to  accompany  these  measurements  with  a  more  developed,  quantitative  microscopi 
theory  that  would  address  how  the  mechanical  energy  of  a  shock  wave,  for  example,  is  converted  m 
chemical  energy  at  the  very  initial  phase  of  the  detonation  process.  Once  this  energy  causes  mi  1 
chemical  processes,  the  mechanism  of  detonation  then  becomes  simultaneously  more  complex  bu 
also  morePin  tune  with  current  theory,  methodologies  and  findings.  However,  it  is  the  microscopic 
nature  of  the  process  at  this  initial  stage  that  is  crucial  since  its  understanding  may  be  expect 
assist  chemists  to  synthesize  new  molecules  of  energetic  materials  that  possess  decreased  sensitivity 

to  deto: shlceitfewS-Sxownthat  different  crystalline  forms  of  an  energetic  material  possess  different 
sensitivities  to  detonation,  the  coUective  behavior  of  the  molecules  in  a  crystal  must  also  be 
understood  for  effective  design  of  new  energetic  materials.  Thus,  the  measurements  and  the  mo 
must  tie  the  molecular  properties  to  those  of  the  aggregate,  the  matrix  or  crystal.  Although  most 
energetic  materials  are  most  commonly  embedded  in  amorphous  solids,  the  crystalline  forms  were 
focuLd  upon  because  they  afford  a  simpler  and  more  quantifiable  system  from  which  a  model  may 
be  built.  Once  formulated,  such  a  model  should  be  capable  of  extension  to  amorphous  composits. 

While  practical  explosives  are  usually  comprised  of  energetic  materials  suspended  m  a  matrix, 
usually  a  polymeric  material,  the  fundamental  microscopic  process  at  the  mitial  phase  of  detonation 
may  be  expected  to  be  essentially  independent  of  the  specific  environment  of  the  molecule.  Under 
this  assumption,  the  model  developed  in  this  research  to  explam  the  microscopic  process  treatst 
surroundings  of  a  given  molecule  in  a  solid  phase  as  being  embedded  in  a  continuous  medium  that 
may  be  represented  by  a  mechanical  susceptibility.  Since  this  mechanical  susceptibility  may  be 
anisotropic,  it  is  represented  by  a  tensor  that  is  itself  related  to  the  elastic  property  ome  n^te^ 

The  experimental  techniques  used  in  the  research  required  large,  well-formed  smgle  ciystals 
of  very  high  optical  quality  for  successful  measurements.  This  proved  to  be  a  major  barrier  to  t  e 
progress  of  the  research  and  to  completion  of  its  goals.  Over  half  of  the  grant  period  was  exhaus e 
in  obtaining  requisite  samples  and  clearing  the  legal  obstacles  for  their  acquisition,  shipping,  s  0  § 
and  study  8  Once  obtained  in  sufficient  quantity,  the  next  difficulty  was  in  developing  a  successfiil 
methodology  for  generating  crystals  of  the  needed  polymorph,  size  ( minimally  2mm  on  an  edge  ), 

optical  quality  and  morphology. 


Stud'es^RDX^u^  ^  prefcrred  method  for  growth  0f  single  crystals  of  RDX  is  by  slow 

evaporation  from  a  saturated  solution  of  acetone.  Using  this  method,  we  were  able  tc .  grow ’smgte 
crystals  but  they  were  small  and  had  fer  too  many  tmy  faces  to  be  useful  in  Raman  and  Brillouin 
spectroscopy.  We  instead  chose  to  grow  RDX  single  crystals  from  a  solution  of  cyclohexanone 
saturated  with  deionized  water.  A  combination  of  gradual  temperature  lowe™8 
resulted  in  large  crystals  with  eight  faces  zonal  to  the  crystallographic  b  axis.  These  crystals  were 
clear,  colorless,  and  visually  defect  free.  Optical  goniometry  revealed  that  two  major  forms,  {00-1 ) 
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and  {100},  that  allowed  easy  study  of  the  different  elements  of  the  Raman  po  amabihty  tensor_  ( 

1  The  previous  polarized  Raman  study  on  RDX  was  performed  by  Rey-Lafon  et  al  m  1 97  • 
They  probed  the  region  from  3200  cm'  to  10  cm' .  There  are  57  possible  internal  vibration  modes 
in  a  RDX  molecule  and  24  Raman-active  lattice  modes.  Rey-Lafon  e/  al  were  able  to  identify 
molecular  modes  and  to  locate  most  of  the  lattice  modes.  Our  work  studied  the  region  from  3300 

cm  to  10cm'.  We  have  located  aU  57  molecular  modes  and  ah  24  of  the  lattice  modes.  Ourmost 

notable  finding  here  is  the  identification  of  a  sharp  fundamental  at  3233  cm'  that  was  previously 
unreportedbyotherresearchers.  Our  instrumentation  allows  us  to  scan  with  a  resolution  greater  than 
1  cnf',  which  aided  in  the  location  of  modes  that  were  apparently  missed  m  the  previous  study. 

Growth  &  Characterization  of  HMX  and  TNAZ 

Single  crystals  of  HMX  were  grown  using  two  methods.  The  first  involved [gradually 

lowering  the  temperature  of  a  saturated  solution  of  HMX  in  acetone.  The  temperature  of  the  solution 
was  controlled  by  thermal  contact  with  a  water  jacket.  The  entire  apparatus  was  thermally  isolated 
to  diminish  any  undesirable  fluctuations  in  temperature.  The  temperature  jncrementaUy 
decreased  at  intervals  of  1  C°  per  day  until  the  appearance  of  small  seed  crystals.  At  this  point  the 
cooling  rate  was  decreased  to  0.1  C°  per  day  until  sufficiently  large  crystals  were  grown.  If  the 
cooling  was  too  rapid,  the  crystals  would  grow  poorly  and  the  process  wou  d  need  to  be  repeated 
The  second  method  involved  an  evaporative  technique.  Rather  than  gradual  cooling  the  saturated 
solution  was  maintained  at  a  constant  temperature  while  the  solvent,  acetone,  was  slowly  evaporated. 
This  method  produced  a  well-formed  crystal  within  approximately  two  weeks.  Compared  to  the 
gradual  cooling  method,  these  crystals  were  grown  relatively  quickly,  therefore,  the  evapora  ive 
method  was  predominantly  used.  Visual  examination  showed  crystals  of  good  optical  quaMy. 

After  the  crystals  were  grown,  their  structures  were  determined  by  X-ray  diffraction 
measurements  that  confirmed  a  non-solvated,  (3-HMX  polymorph,  with  space  group  and  lattice 
constants  consistent  with  those  previously  reported.'  None  of  the  other  polymorphs  were  observed 
The  B-polymorph  crystals  grew  with  a  habit  of  faces  zonal  to  the  “a”  crystallographic  axis.  e  zona 
and  capping  faces  were  determined  with  optical  goniometry  and  later  verified  by  X-ray  diffraction. 

The  crystallization  of  TNAZ  has  been  more  difficult.  Since  evaporation  worked  well  for 
HMX,  a  similar  approach  was  adopted  for  TNAZ.  Initially,  CC14  was  used  as  a  solvent.  As  the 
solvent  evaporated,  no  seed  crystals  would  form  and  eventually,  the  solute  would  solidify  m  either 
a  powder  or  a  poorly  formed  crystalline  plate,  neither  of  which,  is  of  any  use  to  this  study.  Therefore, 
other  solvents  were  studied,  including,  methylene  chloride  and  chloroform.  Although  the  use°ft_h^e 
solvents  afforded  similar  results  to  carbon  tetrachloride,  it  was  observed  that  the  point  at  which  the 
solute  would  solidify  was  very  different.  Methylene  chloride  solutions  could  continue  to  evapora  e 
until  a  highly  saturated  solution  of  TNAZ  would  effectively  precipitate.  However,  when  chloroform 
was  used,  large  plate-like  crystals,  albeit  poor,  would  form  much  earlier.  Therefore  varying  mixtures 
of  chloroform  and  methylene  chloride  were  used.  To  date,  the  most  successful  solution  was 
approximately  80%  chloroform  and  20%  methylene  chloride.  Unfortunately,  the  harvested  crystals 

still  had  too  many  defects  for  any  use  in  spectroscopy. 


Studies  of  P-HMX 

Crystal  Optics  .  .  , 

B-HMX  which  is  a  monoclinic  polymorph,2  requires  a  complete  optical  mdicatnx  in  order  to 

account  for  iny  optical  anisotropy.  To  accomplish  this,  near-normal  incident  specular  reflectivity  was 
measured  at  5 14.5  nm  for  three  symmetrically  inequivalent  faces.  This  wavelength  was  used  for  two 
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reasons  First,  since  the  P-HMX  crystals  are  transparent,  the  imaginary  part  of  the  complex  refrac 
index  can  be  neglected.  Thus,  Fresnel’s  Law  for  reflectivity  is  only  a  function  of  the  refractive  index. 
Second,  this  wavelength  is  the  excitation  wavelength  used  in  both  Raman  and  Bri  own  experiments. 
Thus,  considering  any  dispersion  of  the  refractive  index  with  wavelength  is  unnecessary.  By 
measuring  the  reflectivity  of  light  polarized  along  the  principal  directions  of  each  symmetrica  y 
inequivalent  face,  corresponding  refractive  indices  could  be  calculated  using  Fresnel  s  Law.  In  turn, 
this  data  allows  determination  of  the  semi-axes  of  the  optical  mdicatrix,  which  are:  n,  -  1  58  np  - 
1 .67  and  iy =  1 ,80g  and  from  which  the  refractive  index  along  any  direction  rystal  can  be  calculated. 

Raman  Spectroscopy  of  p-HMX  .  , 

Single-crystal,  polarized,  Raman  spectra  of  P-HMX  were  recorded  from  30  to  3000  cm  with 
varying  polarizations  used  to  assign  modal  symmetries.  The  intramolecular  modes  agree  [see  ig. 

1  ]  with  those  previously  reported.3  The  lattice  modes  assignments  ( spectra  shown  m  Fig S- 2)  “giced 
with  those  reported  by  Iqbal,  et  al?  and  not  those  reported  by  Cavagnat,  etal.  and  Brill,  et  al. 

Beside  these  experimental  studies,  we  have  also  calculated  the  Raman  lattice  frequencies.  T 
calculation  involved  an  empirically  fit,  6-exp  Buckingham  potential  which  is  a  Unction  of  twelve 
variables:  six  lattice  constants,  three  translations  and  three  Euler  angles  that  describe  the  relative 
orientation  of  molecular  inertial  axes  to  the  orthogonal  crystallographic  axes.  The  lattice  energy  was 
initially  calculated  from  these  variables  as  well  as  the  atomic  coordinates  and  then  minimized  with 
respect  to  the  initial  crystallographic  parameters.  The  quality  ofthe  chosen  potential  parameters  was 
judged  on  how  well  the  lattice  parameters  minimize  the  energy  compared  to  experimental  parameters^ 
If  any  ofthe  lattice  parameters  changed  by  more  than  approximately  0.5  A,  the  potential  was  modified 

to  more  accurately  reproduce  experimental  data.  .  , ,  ,  .  ,  . 

Once  the  atomic  coordinates  are  in  the  correct  frame  of  reference,  a  suitable  lattice  needs  to 

be  constructed.  Generally,  a  lattice  consisting  of  150  molecules  will  provide  an  adequate 
representation  ofthe  crystal.  The  lattice  used  for  this  study  was  a  4x4x4  unit  cell  cube.  The  imtial 
potential  parameters  chosen  for  this  system  were  those  of  three  groups:  Sorescu  et  al.  Smith,  etal 
and  Gavezzotti.8  Although  the  first  two  groups’  potentials  focused  primarily  on  mtramine-type 
molecular  crystals,  the  more  general  parameters  reported  by  Gavezzotti  gave  a  minimized  lattice 
energy  of  33  kcal/mol  which  was  the  closest  to  the  experimentally  determined  heat  of  sublimation  of 
41  kcal/mol.9  Although  the  Raman  frequencies  calculated  using  Gavezzotti’s  parameters  did  not 
compare  particularly  well  with  experiment,  his  potential  provides  a  starting  pomt. 

The  potential  parameters  were  modified  through  an  iterative  trial  and  error  process.  This  process 
was  continued  until  the  lattice  constants,  corresponding  to  a  minimized  energy,  were  comparable  o 
experimental  ones.  The  minimized  lattice  energy  was  calculated  to  be  33 .7  kcaFmol  and  the  calculated 
frequencies,  shown  compared  well  with  those  found  experimentally  by  Iqbal  et  a/,  and  us. 

The  results  obtained  over  the  period  ofthe  grant  indicate  that  the  direction  ofthe  research  was 

pursuing  would  lead  to  further  refinement  and  correction  of  data  already  in  the  literature.  Furt  er, 

the  main  purpose  in  obtaining  the  Raman  scattering  data  was  for  its  necessary  combination  with  the 
r  r  _  .  _ xx _ * _ Q+rflin  Gruneisen 


parameters  tnat  are  crmcai  iu  uic  — — - — ■ »  .  t  ~  * 

the  initial  phase  of  mechanical  detonation  of  energetic  materials.  Unfortunately,  the  period  ofthe 

grant  expired  before  the  piezomodulated  Raman  scattering  work  could  be  initiated. 
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HMX  Bg  Symmetry  Optical  Modes 
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Figure  1:  Raman  scattering  of  Bg  symmetry  optical  (molecular)  modes  of /3-HMX. 


RELATIVE  SCATTERING  INTENSITY 


Raman  Scattering  of  Lattice  Modes  of  HMX 


Figure  2:  Raman  scattering  of  lattice  modes  of  0-HMX  for  four  different  scattering  configurations. 


^BriUouMscatterin^measurements  were  made  for  directions  associated  with aU naturally  occimmg 
pairs  of  symmetrically  inequivalent  zonal  faces.  The  low  symmetry  of  most  orgamc  crystals  affords 
only  a  mfoimal  number  of  elastic  constants  to  be  directly  calc, ulated  from 1  aC0lJst^1^S-  J°  j 
monoclinic  system,  the  longitudinal  velocity  for  the  [010]  phonon  is  directly  related  t0  the  C  ^elast^ 
constant.  However,  difficulty  remains  in  the  calculation  of  the  remaining  12  elastic ; 
solution  involves  an  over-determined  set  of  phonons  for  solvmg  the  equations  of  the  Christottel 
determinant  The  experimental  data  for  each  scattering  geometry  mclude  the  direction  cosines  of  the 
phonon  wavevector  the  refractive  index  of  the  incident  and  scattering  waveyector  and  the  measured 
Brillouin  shifts,  which  are  used  to  calculate  the  phonon’s  velocity  m  the  k  -  0  approximation.  Tins 
data  is  then  input  into  a  minimization  program.  An  initial  stiffness  matrix  is  then  incrementally 
minimized  until  the  calculated  values  best  reproduce  the  experimental  ones. 

However,  for  the  scattering  experiments  we  performed,  the  determinant  of  the  calculated  elastic 
constant  matrix  was  not  positive-definite  which  indicates  a  nonphysical  system.  Smce  the  scattering 
eeometriesu^d  involved  pairs  of  faces  zonal  to  the  a  ciystallographic  axis,  scattering  was 
constrained  to  the  be*  plane.  Based  on  the  poor  convergence  of  the  stiffiiess  matrix,  the  elastic 
anisotropy  of  HMX  cannot  be  folly  described  using  phonons  of  a  single  plane.  Therefore,  adaptations 
to  the  Brillouin  spectrometer  were  made  that  facilitated  the  measurement  of  non-zonal  capping  faces 
Besides  simply  increasing  the  number  of  scattering  geometries,  using  p honons of 
permits  accounting  for  the  anisotropy  of  the  elastic  constant  matrix.  However,  the  period  of  the  grant 
expired  before  the  spectra  on  the  modified  spectrometer  could  be  obtained.  , 

F  Although  the  stiffiiess  matrix  could  not  be  calculated  directly  from  experimental  data  an 
approximate  set  of  elastic  constants  can  be  calculated  using  theoretical  scattering  data.  The .direction 
cosines  of  the  phonon  can  be  accurately  calculated  for  any  scattering  geometry.  Smce  the  HMX 
optical  indicatrix  was  measured,  the  refractive  index  along  any  incident  or  scattering  direction  cam foe 
calculated.  This  leaves  the  Brillouin  shifts  that  can  be  well  approximated  by  comparison  to  the  large 
number  of  actual  spectra  taken  from  phonon  scattering  in  the  be*  plane  From  these  Brillouin  shift  s 
we  have  been  able  to  calculate  the  approximate  stiffiiess  matrix  (Fig.  3)  which  can  also  be  used  to 
calculate  such  properties  as  the  bulk  compressibility  and  Young’s  modulus. 
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Theory  of  Initiation  of  Detonation 

Over  the  period  of  the  grant,  collaboration  with  Prof.  Tadeusz  Luty  of  the  Technical 
University  of  Wroclaw  was  actively  focused  on  formulation  of  a  quantitative  microscopic  mo  e  or 
detonation.  Given  the  breadth  of  the  subject,  the  focus  was  specifically  on  how  mechanical  energy 
is  directly  converted  into  chemical  energy  at  the  very  initial  stage  of  detonation  where  the  shock  wave 
traverses  the  unreacted  and  previously  unperturbed  molecules  in  the  crystal.  A  microscopic  mo  e 
dealing  with  this  phase  of  detonation  may  be  expected  to  be  useful  to  synthetic  and  materia  s 

chemists  in  the  design,  synthesis  and  production  of  new  energetic  materials. 

The  model  that  was  developed  is  based  on  the  idea  of  Gilman  that  a  molecule  does  not  ha 
time  to  mechanically  respond  to  the  shock  wave  through  its  vibratonal  degrees  of  freedom  but  rat  er 
must  take  up  the  available  energy  by  closing  the  highest  occupied  molecular  orbital-lowest  occupied 
molecular  orbital  (HOMO-LUMO)  gap.  However,  it  is  less  restrictive  and  does  not  require  complete 
closure  of  the  gap  and  subsequent  metallization  to  explain  how  the  mechanical  energy  is  coverte  o 
molecular  excitation  and  reaction.  The  theory  shows  that  the  deformation  potential  of  the j  molecule^ 
<E>def,  is  the  crucial  physical  property  and  that  this  determines  the  degree  to  which  the  HOMO-LUMO 

gap  must  decrease  to  lead  to  an  initial  chemical  response  to  the  shock  wave. 

Our  theory  for  mechanochemical  induced  processes  was  published  m  the  J.  Chem.  Phys. 
The  treatment  shows  that  the  initial  response  of  a  molecule  to  a  shock  wave  requires  the  molecu  e 
to  distort  to  a  geometry  that  is  essentially  the  “ground  state”  of  the  molecule  suffering  the  extremely 
high  effective  pressure  of  the  shock  wave.  A  new  insight  in  this  treatment  is  the  representation  of  the 
distortion  of  the  molecule  by  a  linear  combination  of  normal  mode  coordinates.  This  compressed 
ground  state”  arises  from  mixing  of  the  ambient  ground-state  potential  of  the  molecule  in  the  crys  a 
with  a  dissociative  excited  state  potential.  The  combination  of  these  two  diabats  to  form  the  high- 
pressure  ground  state”  is  represented  in  Figure  4. 


(M+M) 


LUMO-f 

HOMO+ 


-LUMO 
4+HOMO 


<Q>C  <Q>„ 

Figure  4:  Diabats  for  mixing  of  states  for  a  molecule  under  significant  stress.  Qc  is 
the  normal  coordinate  representing  some  critical  deformation  while  Qd  is  the  value 
of  that  coordinate  for  the  deformation  mixing.  LUMO  and  HOMO  are  the  Lowest 
Unoccupied  Molecular  Orbital  and  Highest  Occupied  Molecular  Orbital,  respectively. 
I  is  the  ionization  potential  and  A  is  the  electron  affinity.  M  represents  the  molecu  e 
in  various  valence  states.  4>def  is  the  deformation  potential  and  t  is  the  energy  of  the 


gap  between  the  new  states  arising  from  mixing, 
product  and  reactant  molecular  species. 


E  and  ER  are  the  energies  of  the 
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The  height  of  the  barrier,  3>def,  is  related  to  how  easily  the  molecule  can  be  deformed  and  thus,  is 
obviously  specific  to  a  given  molecule.  If  there  is  excess  mechanical  energy,  then  the  molecule  can 
be  forced  over  the  barrier  to  produce  products  that  will  be  essentially  thoseassociatedwithth 
mixing  dissociative  state  but  there  will  be  products  that  also  relate  to  thermal  decompositon  of  he 
bound  ground  state.  Thus,  the  mechanical  energy  is  used  to  deform  a  molecule  to  yield  the  dlustra ted 
potential.  If  there  is  excess  mechanical  energy  beyond  that  required  for  this  deformation,  the 
molecule  can  traverse  the  potential  barrier  to  produce  products  consistent  with  those  of  the  wo 
diabats.  The  model  thus  explains  why  photochemical  products  of  the  energetic  materials  are  observed 

as  understand  the  process  is  by  analogy  to  the  reaction  rate  theory  of  Marcus. 

In  the  illustration  [Fig.  5]  the  mixing  of  the  “reacting”  states,  in  this  case  the  ground  state  and  the 
excited  state  that  is  accessed  because  of  compression  driving  the  molecule  up  the 
its  ground  state  potential,  is  shown  on  the  left  in  common  representation  but  onthe  left  the  plot  is  that 
of  the  resulting  high  stress  “ground  state”  represented  in  a  reaction  coordinate  plot. 


Figure  5  Representation  of  compressive  excitation  of  ground  state  to  intersection  with  higher  energy 
state  where  mixing  may  occur  (left)  and  showing  reactive  nature  of  the  resulting  high  pressure  ground 
state"  as  a  reaction  coordinate  diagram  (right). 


PThe  role  of  “hot  spots,”  regions  of  the  crystal  of  the  energetic  material  that  appear  to  be 
involved  “prematurely”  in  the  initial  detonation  process,  has  been  of  interest  for  some  time  The  idea 
that  such  “hot  spots”  are  associated  with  lattice  defects  has  long  standing.  This  model  gives 

quantitative  expression  to  this  qualitative  idea.  ,  •  „ 

Because  defects  are  invariably  associated  with  a  strain  field,  molecules  in  the  region  of  the 
defect  are  already  under  some  tension,  perhaps  enough  to  distort  them  to  the  stress-induced  state  of 
Fig  5  (right)  Such  molecules  are  already  “on  their  way”  to  passing  over  the  barrier  and  subsequent 
formation  of  detonation  products.  The  height  of  this  barrier  is  the  deformation  potential  and  since 
the  strain  field  is  already  contributing  to  the  distortion  of  the  molecules  in  its  region,  such  molecules 
are  already  possessed  of  the  energy  needed  for  reaction.  Alternatively,  even  if  the  molecule  is  no 
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initially  in  the  “compressed  ground  state”  it  is  already  in  a  sufficently  strained  region  that  it  does  no 
require  as  much  energy  as  a  molecule  in  a  strain-free  region  requires  to  satisfy  the  need  of  sufficien 
energy  to  surmount  the  4>dd- barrier.  Such  molecules  will  thus  be  the  first  to  decompose  Either  way, 
the  model  accounts  for  “hot  spots”  and  substantiates  the  view  they  are  associated  with  defects. 

“Inverse”  Jjihn-Tcllcr  Effect  i  j 

A  further  result  of  the  mechanochemical  theory  is  that  symmetry  arguments  based  on  the 

equations  show  that  high-symmetry  molecules  are  necessarily  more  unstable  to  mechanical  stress. 
This  may  be  thought  of  as  the  molecule  reacting  to  the  imposed  stress  but  unable  to  accommodate 
the  stress  by  deformation  and  thus  only  being  able  to  achieve  relief  by  bond  scission.  This  has  been 
seen  in  a  full  ab  initio  calculation'2  that  shows  RDX  breaks  only  „neN-N  breaks  m  the  mitml  process. 
We  were  able  to  confirm  this  by  a  DFT  calculation  based  on  our  modef  However,  the  operant 
principle  is  that  of  symmetry  and  it  has  a  direct  message  for  chemists  regarding  the  role  of  symme  ry 

in  the  design  of  energetic  materials  for  explosives. 

Although  the  subsequent  fate  of  the  energy,  formation  of  the  detonation  wave,  etc.  was  not 

addressed  in  this  theoretical  work,  the  aspects  that  may  be  argued  to  be  most  useful  to  the  sjmthetic 
chemist  trying  to  design  a  new  explosive  are  delineated.  This  can  only  be  achieved  from  a 
microscopic  understanding  of  all  of  the  processes  involved  m  detonation.  This  particular  model 
produces  a  “zeroth-order”  microscopic  picture  that  can  be  used  in  rational  design  of  energe 

materials. 
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